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ON THE SMALL DIVISORS IN THE LUNAR THEORY* 

BY 

ERNEST WILLIAM BEOWN 

Introduction. 

1. One of the greatest difficulties which occur in the actual calculations 
necessary to obtain the coefficients of the periodic terms which express the co- 
ordinates of the Moon in terms of the time, is the presence of small divisors. 
At each stage of the approximations either one or two such divisors enter in 
each set of terms,f and the entrance of a small divisor means, in general, a loss 
of accuracy in the corresponding coefficient. To illustrate the manner in which 
these divisors enter, it is sufficient to refer to the general forms which the equa- 
tions of motion take. For example, consider the equations 

X + x = '^a. cos (2^ -|- a) i , 

< ' (f = 0, d=l,±2, •■•), 

y = ^6. cos (2i-f a)t 

i 

where a^ ,h., a are constants, a being nearly unity or nearly zero. 
The particular integrals are 

The small divisors occur for 

i= 0, —1 In X when a = 1 nearly, 

■1 = in y when a = nearly. 

The corresponding coefficients are increased in size. If a numerical develop- 
ment is being made, the number of places of decimals to which the results are 
accurate will be less in the integrals than in the differential equations. If a 
literal development is under consideration such that a., h., a are series in posi- 
tive integral powers of a small quantity m , and if a or 1 — a be divisible by 



* Presented to the Society February 22, 1903. Received for publication February 4, 1902. 
t A. 'set of terms ' is defined in Article 10 below. 
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m , the particular coefficients referred to will be known to one less power of m 
than the corresponding coefficients in the differential equations. 

2. Now it is well known that, for all inequalities depending on the attractions 
of the sun and earth alone,* whatever small divisors may occur, the results 
never contain negative powers of m , where m is the ratio of the mean motions 
of the sun and moon.f Hence, if a small divisor contains the factor m,'', the 
corresponding coefficients a^, h^ must also contain this factor. 

The coefficients of the periodic terms in each approximation depend on the 
results of previous approximations and there is, in general, a loss of accuracy at 
each stage, which arises in the following way. Suppose that a coefficient in the 
differential equations appears in the form 

(a) {\ + \m + \m^ + ■■■ + h^m") - {\ + h^m + h'^rn'+ ■■■ + h^m"), 

where the h^, li\ are rational numerical fractions, i. e., it appears as the differ- 
ence of two expressions which are known to m"". If the divisor be m,, the re- 
sult is only known to irf~^. 

If, however, the calculations could have been so arranged that the coefficient 
appeared in the form 

(/3) {\m-\-\m'+---+h^m'){h[m + h'^m'+ ... + 6;™"), 

the result after division by m, would still be true to m'',\ that is, there would 
have been no loss of accuracy. 

3. If the numerical value of m be used from the start, similar results hold ; 
instead of considering powers of m , we consider the number of places of decimals. 
The question of accuracy is a little more complicated since the coefficients of the 
higher powers of m are frequently large numbers, that is, the series arranged in 
powers of m converge slowly in many cases. But the gain in accuracy is great 
since m only occurs in the original differential equations as a factor to the first 
and second powers, so that we exchange a slowly converging (or perhaps diverg- 
ing) series in powers of m for a series ^ 6^10~' where h^ is a positive integer less 
than 10. 

In the lunar theory m is about -^^ and the small divisors are generally infinite 
series of positive powers of m , beginning with the first power. The loss of 
accuracy caused by such a divisor at each stage of the approximations is gener- 

*The ratio of the mass of the moon to that of the sun is supposed to be neglected. 

fin this theory m is the ratio »'/(" — "')• where n', n are the mean motions. The state- 
ment is true for either m or m. 

JThe phrases "true to m*" or " accurate to m* " are taken throughout to mean that all coeffi- 
cients up to and inclusive o/that of m'' are accurately known. A similar meaning is attached to 
the phrase "accurate to p places of decimals." 
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ally balanced by the presence of a new small factor independent of m and smaller 
than m ; this factor is one of the four quantities : e = the eccentricity of the moon's 
orbit, e' = that of the sun's orbit, k = the inclination of the moon's orbit, a = 
the ratio of the mean distances of the moon and sun. The numerical accuracy 
would be the same at any stage of the approximations as it is at the start were 
it not for two circumstances. First, at every stage of the approximations after 
the second (that is, when the order with reference to e, e, k, a is greater than 
the first) there are one or more terms having as divisors series which begin 
with m^ ; second, the numerical multiplier of the lowest power of m present in 
a coefficient may be large, and these multipliers tend, in general, to increase 
with the order of the term. These facts require the earlier approximations to 
be taken to a higher degree of accuracy than would otherwise have been neces- 
sary. There are no divisors beginning with m' within the practical limits set by 
observation.* The loss of accuracy due to a large numerical multiplier of the 
lowest power of m is judged by the ratio of the numerical value of the coefficient 
to Xm' where X is the characteristic f and m** the lowest power of m present. 

When the literal value of m is used all the multipliers of the different powers 
of m are rational numerical fractions. If the numerical value of m be used from 
the outset, m being a quantity obtained from observation to a certain number of 
places of decimals, there is a loss of accuracy due to the use of a finite number 
of places of decimals in long and complicated calculations. I have examined 
elsewhere :j: the maximum loss which can be produced from this cause and have 
shown that when the number of operations (addition, multiplication, etc.) is 
less than 63 the loss is not more than two places of decimals. As a matter 
of fact, however, the actual loss will never be more than this, however ex- 
tended be the calculations, since positive and negative errors tend to balance 
one another. 

4. The main object of this paper is to show how the method § now being used 
for the calculation of the inequalities produced by the attraction of the sun can 
be modified so that tfie total loss of accuracy, as far as the small quantity m 
is concerned, does not exceed q — 2 porcers of m when the -order of the char- 
acteristic is q, with the exception of certain coefficients of the second and third 
orders. 

This theorem is proved for a literal expansion in positive powers of 
m, e, e', k, a. From what has been said above it is not difficult to make an 

* The largest coefficient with such a divisor has the factor ek^c' m at least. In Delaunay's 
notation this term lias the argument I -\- 21-' — 41) — At', and the loss of accuracy for it is only 
one power of m. Cf. Art. 14, etc. 

fThe characteristic is defined in Art. 10. 

t Hdverford CnVege Studies, no. U (1892). 

§ T. M. M., part I. See footnote 1 on page 164. 
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estimate of the degree of numerical accuracy required at the outset for a given 
degree of numerical accuracy in all coefficients. 

First, e, e', k, a are all less than -^^, so that a characteristic of order q con- 
tains a factor 10~* at least. Second, large numerical multipliers cause a loss of 
two places of decimals at most within the limits set by observation. Third, an 
estimated loss of two places will occur owing to the use of decimals. Hence, 
within the limits set by observation, this maximum loss of accuracy in the terms 
oforderyis ^-2-^ + 2 + 2 = 2 

The two quantities which require to be most accurately known for comparison 
with observation are the mean motions of the perigee and node and it is useful 
to have the ratio of these to the mean motion of the moon accurately to 9 
places of decimals. It is not necessary to know the numerical coefficients of 
the periodic terms to more than 8 places. The intermediate orbit — Hill's 
variation curve * — used in the theory is the first approximation. In finding the 
principal parts of the motions of the node and perigee there is a gain of one 
power of m,j- and an estimated loss of two places of decimals owing to ex- 
tended calculations over the results in the variation orbit — a net loss of one 
place. Hence, for these mean motions we require the variation orbit (which 
contains no small divisors) to 9 -|- 1 places, and, from the statement above, for the 
the coefficients of the periodic terms it should be known to 8 -f 2 places. 
Hence, for all purposes it is sufficient to know the coefficients in the intermed- 
iate orbit to 10 places of decimals. 

6. In order to clear the way for the proof of the theorem enunciated in Art. 
4, a general formula is given in Sect. I for the gain or loss of accuracy in such 
expressions as actually occur in the theory. In Sect. II will be found the differ- 
ential equations in various forms and an investigation of the expressions which 
are affected by small divisors. Sect. Ill contains certain transformations which 
assist in preventing a loss of accuracy. In Sect. IV the detailed discussion 
necessary to prove the theorem for the first four orders is given with a general 
indication of the proof for higher orders. 

It is quite possible that transformations might be found such that no loss of 
accuracy, as far as powers of m are concerned, occurs. Such transformations 
would doubtless be somewhat complicated and difficult to arrange for computa- 
tion. Moreover, from the result at the end of Art. 4, they are not practically 
necessary and it does not seem useful to extend analytical transformations, already 
somewhat long, beyond a point where theoretical or practical interest ceases. 

*G. W. Hill, Researches in the Lunar Theury : American Journal of Mathematics, 
vol. 1 (1878). 

fG. W. Hill, Literal expression for the motion of the Moon's Perigee: Annals of Math- 
ematics, vol. 9 (1894). The same result is obviously true for the motion of the node. 
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Section I. A criterion of accxjract. 
6. Consider the expression 

ZS1S2 ■ ■ • Sp_ _ n >s'j 

where S^, T^ are power series in m, each of them known up to some definite 
power which may differ in the different series, the coefficients in the series be- 
ing rational numerical fractions independent of m, and where X) denotes a sum 
of such expressions. It is desired to find the highest power of m whose coeffi- 
cient is accurate when the expression is expanded in powers of m . 

As there is no need to consider the numerical coefficients in the power 
series, here and in the following sections (a , /3) will be used to denote a power 
series in which the first power of m whose coefficient is not zero is m" and in 
which the greatest known power of m is m^. In accordance with this notation, 
put 

It is desired to find B . When B is known for each such expression, the re- 
quired number is the least of the numbers B present in ^ . The second num- 
ber in any parenthesis will be called the degree of accuracy. 

Suppose that any one of the series is known to any power of m , i. e., that 
the coefficients proceed according to a known law, or that its accuracy is greater 
than any we need to consider ; we can signalize the fact by putting oo for the 
second number in the brackeb. Then 

(a , OD ) = m" (0 , OD ) , 

and as (0 , od ) does not affect the accuracy of the result, we can collect all such 
terms in the numerator and denominator and consider the expression 

First, noting the lowest power of m occurring in each series, we have 

A = ^a-Y.l + 1c. 

Second, the number of terms in (^A , B) cannot be greater than the number 
of terms in any one of the component series. Hence 

B — A^ the least of /8; — a^ , h^ — <y^. 

These results determine A, B. 



164 
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7. In the applications to be made the following limitations may be imposed : 

(i) For all factors of the denominator in which 7^ = , S. — y.is equal to or 
greater than any of the other B. — 7^ or than any of the 0. — a^; 

(ii) The factors of the denominator in which 7^ 4= , are at most two in num- 
ber and then of the form (7, BJ^; 

(iii) a., 7 are positive; 

(iv) 7 < 3 . 

Thus the expression reduces to 



and in this 
(7) 



m*(aj,/3J...(a^,/SJ-(7,S)*=(^,5) (2 = 0,1,2), 

A= k — q'Y+ J^a, 
B — A = the least of the numbers 0^ — a., S — 7 . 



(v) The numbers a., /S^ , A, B will denote the numbers set down or any re- 
spectively greater numbers, but 7, B always denote the definite numbers set 
down. 

Section III. The differential equations and their solution by 

APPROXIMATION. SmALL DIVISORS. 

8. The differential equations. The notation adopted in the series of mem- 
oirs on the motion of the moon* will be in general used here. The equations 
of motions are f 



(1) 



a. 



h. 



c. 



<JS 

D's - 2m Z>s = f^, 

au 



D^z 






where x , y, z are the coordinates of the moon's center and 



u = x + yV—'i-, s = x — yV—'^, r^ = x'' + y^ + z^ = us + z^, 

n 



D = - 



V-\ d 



m = 



n — n dV 

n, n being the observed mean motions of the sun and moon. Also 

* Ernest W. Brown, Theory of the Motion of the Moon, Memoirs of the Eoyal Astro- 
nomical Society, part I (1897) in vol. 53, pp. 39-116 ; part II (1899), pp. 163-202 ; part III 
(1900) in vol. 54, pp. 1-63. These will be referred to under the abbreviation "T. M. M., pts. 
I, II, III." 

fT. M. M., part I, section (ii). 
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/ =/-! +./; +/3 + • • • = Ep/, (p = -l,2,3, 4, •■■), 



(Id) 



/j = — |m^(M + sf + m^a^, f2 = e'm^ (function of m, s, ^), 



in which y is a homogeneous function of m , s , s of degree p containing only- 
even powers of z , and always multiplied by the factor m^ except when p = —1. 
Denoting by «„, a' the mean distances of the moon and sun used in the theory, 

«=(l+2m+ ...)al, 

and^ contains the factor {aja')''~^ = a'"'^ when p > 2; y also contains t ex- 
plicitly when p = 2. 

9. From the equations (1) the following are easily derived :* 

a. uDs — sDu — 2mus = Z)""' I m ^ s-^ j , 

B. DuDs+ {Dzf ^ C+f - Z>-> (Z>/) , 

C. D\us + »' ) - 2Du Ds-I {Bzf - 2m (uBs - sBu) = X) pf^, 
d. B" {us + z^)-Du-Ds- {Dzf - 2m {uDs - sDu) 



(2) 



e. DiuDz — zDu) - 2mzDu = lu %- z |^, 

^ oz ds 



f. D\us + »2) -DuDs- {Dzf + uDs - sJDu - 2m(2m + 1)ms 

= a+i:(p + l)/,-t-i)-{(2m+l)(r.^-3^)-i>'/}, 
^. Z>2(ms+»2-)_2),^.2>s_ (Z)»)2+ 1 (^<Z>s-sZ>M)(l-m)-m(l + 3m)Ms 

where Z)y denotes the derivative of f with respect to i/ — l(ra— n')^ only in 
so far as t is contained explicitly iny, D~^ denotes the operation inverse to D, 
i, e., an integration with respect to j/ — l{n — n' )t, and C is a constant of 
integration. The equations {2B), (2C), distinguished by capital letters alone 
contain y_j . If the parallax of the sun be neglected f=f_^ +f^ and the equa- 



*T. M. M., pt. I, pp. 53, 54, contain equations (2a) — (2e) ; (2/), (2y) are obvious conse- 
quences of these. 
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tlons with small letters are all of the second degree with respect to u,s, z ex- 
cept for the constant C; if the eccentricity of the sun be neglected, D'f= ; 
if both these constants be neglected/" =y_, + f\ and the equations give all 
terms depending only on m and the eccentricity and inclination of the moon's 
orbit. Of these seven equations three only are necessary to replace equations 
(2) and one of the three must be (2e) . 

10. The values of m, s, z derived from these equations are known to be for- 
mally expansible in cosines and sines of multiples of the angles (n — n'^t -\- ^|,, 
n't -f i,, c(ra — n')i + t^^ ^{n — n'^t -\- ^3, with coefficients expanded in ascend- 
ing powers of m, e, e', k, a, where c, g are constants to be determined. If we 
put 

the periodic terms can be expressed in positive and negative powers of h^ ^, h^ ^'", 
h^^°, SjS'S where &^= e''^"''- = const. Since we are only concerned here with 
the determination of the coefficients of the periodic terms and since m, c, g are 
supposed to be incommensurable with unity and finally, since the method used 
is that of indeterminate coefficients, the constants h^ may be omitted throughout 
and then all periodic terms are expressible in positive and negative powers of 

The general expressions for u ^~^ , z l/— 1 are then sums of terms of the form 

a„ e*i e'*2 k*' a** (X. §"^+2' + \'_. };-'-ii) , 
where t = (^i\ — i^)c + (2%.^ — i^)m + {2i'^ — i^)g; 

V ^-^ ^' ■'■' ^' ■■■' ^^--v 

2^ = , =fc 1 , ± 2 , • • • ; 2i — ^^ = even no. ; 

in w?~S ths values of i^ are even, 

in z \/ — 1 , the values of i.^ are odd and X^^ = — \; 

Q^ is the linear unit; X^, X^^ are positive power series in m for which the nota- 
tion * 

f f f 

with Xj obtained by interchanging e with e', rj with rj', k with k' , is used except 
when t., %'. are all zero when X, = a,. 

A set of terms contains all the terms for which j , i' are the same, i. e., 
the sum when i only receives all values. The characteristic of the set is 



*T. M. M., pt. I, p. 60. The form in which the notation is presented here has been slightly 
changed but its essential character is unaltered. The coefiBcients for terms of the first order are 
(«)i, {'' )i, • • • , («)i ; these are written f^ , £^ , • • • , Oj, for the salse of brevity. 
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e''e"''k'«a'<, and the order of the characteristic or simply the order of the set 
is ij + i^ + *3 + ^4 • The argument of any term is the index of f in that term. 
The calculations show that 

c=l + m — fm'' , gr = 1 + m + |m^+ •■•, 

where the remaining parts consist of higher powers of m and terms involving 
powers of m and powers and products of e?, e , k^, a^. For our purpose, the 
first two terms in c, g are alone of importance; the others all contain the factor 
m^ (see Art. 3). 

11. The method of approximation consists in successively determining the 
terms of orders 0, 1 , 2, • • •. A set of terms of any given order can be found 
independently of all other sets of the same or higher order, but the coefficients 
in any one set depend on one another in the forms which the equations take 
here,* so that it is necessary to consider all the terms of a set together. Let 

w = w„ + Wj + 1*2 + • • • > « = 2i + 2, + • • • , 

where u,^, z,^ contain all terms of order X, and suppose that we wish to deter- 
mine terms of order X (X 4= 0) in m or is , those of lower orders having been 
obtained. 

Let Q(u, s, z) be any one of the functions which we shall use, and let [QJ^ 
denote the terms of order X in § after expansion has taken place. Then 



Put 



[/ A— 1 A-1 A-1 \-| 

so that, after expansion by Taylor's theorem. 

The known terms are thus separated out from the unknown terms and all the 
equations will be linear with respect to the unknowns. All our functions are 
such that Q contains either only even or only odd powers of a; in the former 
case the fourth term, in the latter ease the first three terms are absent from [ § ]\ . 
When this notation is applied to the equations (2) the remarks of Art. 8 show 

th at [D-'( Df)-],^lD-\Df)]:, [/Ja=[/JI for p>2. 

*In a paper On the solution of a pair of simultaneous differential equations, etc. (Memoirs 
of the Cambridge Philosophical Society, vol. 18 (Stokes Memorial, 1900), pp. 94-106), 
I hare shown how any coefiBoients in the set may he determined independently of all other 
coefficients of the same set when the terms of lower orders have been found. 
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The symbol P will be used as an abbreviation of the phrase " a positive 
ascending power series in m beginning with the gth power of m." The letter 
P (without a suffix) stands for a similar abbreviation with the addition " 5'= 1 ." 

12. The first approximation consists of the terms of zero order and is given by 

«„?"' = Ej«jr' (t=o, ±i,±2, •••), 

in which a^ is a certain length which we shall here put equal to unity, and 

/„=-2-4m + P„ uf(^ = l + ^m + P^, s„ |^ = 1 + 2m + P,. 

The set of terms with a characteristic of order X and arguments ± (t + 2i) 
are, for te, given by , , , ,„. , „. 

where the characteristic is omitted. Substitute 

A A-l 

in equations (2) except (2e) and equate to zero the coefficient* of ^'■+^*. The 
various terms in the coefficient may be divided into three classes. The first class 
consists of the terms containing a^X^, %X'_^ , or V, \'_^ since a„ = 1 ; the second 
class consists of terms of the form \P„i ^'jPq{l~^t j = ^ 1 ^ 1 ) =t 2 , • • •) , 
since a. = P„.r, the third class consists of known terms which do not involve 
X., X'.. The terms in the first class are those of principal importance in the 
determination of X^, Xl., since on them depends the existence or non-existence 
of small divisors. 
For brevity writef 

X,X',T instead of X^, X^^, t -f 2i; 

then the terms of principal importance in equations (2) are respectively 

a. _t(X-X')-2(1 + m)(X + X'), 

B. _T(X-X')-2(l + m)(X + X'), 

a 2T(H-m)(X-X') + (T^+3 + 6m)(X + V), 

d. T(l + 2m)(X-X') + (T2+l+4m)(X + X'), 

/. 0(X-X') + (T2-l-2m)(X+X'), 

g. i(l+m)T(X-X') + T^(X + V). 



(3) 



*The coefficients of f=^('^+2*) are equal in the equations. 

t The fact that the symbol 1 has been previonaly used for the characteristic does not cause any 
confusion, and these abbreviations are of considerable advantage. 
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These results show that if any pair of the equations (2) , excluding (2e) and 
the pair (2a) , (25) , be used, the determinant of a pair of equations for finding 
X , \' is , except as to a constant factor which is not small, 

T(T2-l-2m), 
and this is equal to jP if either 

'r=P,^ ±T=l + m+P^ (2 = lor2). 

I have shown that the more accurate statement* of the case is that smaU 
divisors equal to P occur if, and only if, 

but the above statement is sufficient since c = 1 + m + P^ and since we do not 
consider values of q greater than 2. 

13. Three theorems will now be proved, two of which are necessary for the 
discussion which follows. 

Theorem I. If±r=.\ + XQ.+ P, then 3X + X' can he determined with- 
out loss of accuracy due to the small divisor. 

For by (3a), (Bd), the equations (2a), (2d) determine 

- (1 + m)(3X + X') — (X — X') (t — 1 - m) , 

(1 + 2m) (3X + X') + (1 + 2m) (X - X') (t - 1 - m) + P^' 

respectively. Hence multiplying (2a) by 1 + 2m and subtracting from (2d) , 
we shall have an equation to determine (1 — m) (3X + X'). AU the other terms 
containing X, X' in these equations have the factor P^ at least, and therefore, 
since q'^2, there is no loss of accuracy in the determination of PjX, PjX' and 
therefore none in that of 3X + X'. 

Theorem II. If t = P , then X + X' can he determined without loss of 
accuracy due to the small divisor. 

For (3/), (2f) show that the first approximation to X + X' has this property 
and (3P) , (2P) show that then X — X' has the divisor t . But the other terms 
in (2f) containing X, X' have the factor P^ at least, and there is no loss of ac- 
curacy in PjX or PjX '. 



* Investigations in the Lunar Theory, American Journal of Mathematics, vol. 17 (1895), 
p. 356. This statement involves the divisor r^ = P^ when r = P^ . Moreover, the right hand 
member oE equation (2a) involves the divisor r. But whenever ?., V are contained in the right 
hand member of (2a) they occur with the factor m^a±i -;- t = P4_j which is equal to P^ in the 
worst case (g = 2) and we are neglecting P^ in the terms of principal importance. It might still 
be thought that when g = 2 , the terms involving m' might disappear in the determinant, so as 
to make the divisor equal to P^. That this is not the case I shall show in a future paper. 
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Theorem III. All the coefficients in l/r for which t = Pj are equal to 
(2 , /3) and those for which r = P^ are equal to (1 , /3). 

If r = P^, (2a) shows that uDs — sDu = (1 , /8) ; (2d) then shows that 
I)u-Bs+{I)zf={2,^); {2B) then shows that /= (2, /3). The result 
follows since 

/ = /_! + (2, ^) = (1 + 2m)/r + (2, ^). 

If T= P^, (2d) shows that Du ■ Ds + {Dzf = (1 , /8) and thence, as before, 
from(25),l/V=(l,y8). 

These three theorems show that the functions which may be found without 
loss of accuracy are 

3X + X', (t'-c')(X + X') for±T=c + P^; 

(9= lor 2). 
X + X', t(X - X') , l/r = P^_^ for t = P^ 

The degree of accuracy here referred to is that to which the known terms in 
the equation have been previously found. 

14. The special portions of the equations whose accuracy is affected hy the 
small divisors. There are four cases to consider 



T 



= P. T=l+m + P„ (g = lor2), 



? 



The case — T=l+m + P is the same as t = 1 + m + P since the change 
— T for T merely alters the sign of X — X'. In all cases C is absent. 

(i) T = Pj . Equation (2g) determines t (X — X') without loss of accuracy. 
Also the right hand member contains the factor m^ -;- t ^ P^ at the worst and 
is therefore divisible by t without loss of accuracy. The only terms in the 
right hand member which are not divisible by t without loss of accuracy are 

(4) [- DuDs- {Bzf + -|- (uBs - sDu)'] ^ 

(4)' = [- i)(M?-')-Z>(sO - J«?->2)(sr)+ isri> («?-') - (i>«r]x; 

for IP (us + s^) gives the factor t^ and the terms other than those just written 
have the factor m . 

The problem here is so to express (4) or (4)' that m or P^ appears as a factor 
explicitly ; on division by t, the degree of accuracy then remains the same pro- 
vided this factor itself is sufficiently accurate. 

(ii) T = Pj. It would here be necessary to transform all terms in (2^) except 
D^(us + a^), 3m^Ms, "^{p + 1)/ which contain factors t^, m^, m^ respectively, 
for /)"' produces a divisor P^, and the final division by t produces another such 
divisor. A special transformation is given in Art. 17 for this case; it will also 
serve for the previous case. 
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(iii) T — c = Pj or Pj . Equation (2/") determines 

(t^ _ 1 _ 2m)(X + V) = P^(X + X') (2 = 1 or 2) 

without loss of accuracy. The right hand member has the factor m^ and is 
therefore divisible by P^ without loss of accuracy. The portion of the left 
hand member which does not contain m^ explicitly is 

lD\us + z')-Du-Ds- {Dzf + uBi - sDu - 2m«s]^, 
or, since D^ produces t^ and t' = 1 + 2m + P , 
(6) l-DuDs + uDs - sDu + us - {Dzf + Z>^ (s^)] ^ 

(5') = - [Z> («?-! ) • D (.?) + {Dzf - D' {z' )] , , 

the terms which are explicitly divisible by t — c = P being omitted. This 
last is the expression to be transformed so as to show the factor P^ explicitly. 

15. The previous results only refer to the terms in m; we have still to consider 
the equation for z. Consider equations (Ic), (2e) and choose out the terms of 
order X. The terms of principal importance for the determination of the co- 
efficient of f ■^ are obtained by substituting 



^,i/-i = x(r^-?-), « = «„, 

and equating to zero the coefficient of ^=^'' in (ic) or of ^*''+i in (3e). Powers 
of m beyond the first being neglected, the equations give, for the terms of prin- 
cipal importance, 

(T2_l_2m)X. 

There ia thus a divisor P when t = 1 + m ± P . The same remark as that 
made in Art. 12 applies here also, but with this difference: that the small di- 
visor approximates to t^ — g^ instead of t^ — c^. This difference, however, does 
not affect our argument, since c, g only differ by P^ and we are not consider- 
ing divisors equal to P3. 

16. I shall now prove the 

Theorem. — All coefficients in z can he determined without loss of accuracy 
due to small divisors in the equation for z. 

The only cases where there may be a loss of accuracy are In those coefficients 
where ± t = 1 -f m -f P^ (§■ = 1 or 2) . 

Consider equation (3e). The only terms which do not contain m^ explicitly are 

(a) [Z> {uDz -zDu) — 2m zDu']^ . 

1°. Suppose T = 1 -I- m -f P . Equate to zero the coefficient of ^'"'■. 
Then (a) becomes, for this coefficient. 
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[_ (m + P^{uDz - zDu) - 2m»Z>M];,, 

= [— P (uDz — zDw) — mP («»)] X 

= [- P^{uPz - zPu) + m(m + P^)uz]^ . 

On dividing this by t^ — 1 — 2m = P + P^ no accuracy is lost since P or 
m^ are present explicitly. 

2°. Suppose T=— 1 — m + P. Equate to zero the coefficient of ^'+^. 
Then a becomes of the same form as before and consequently no accuracy is lost 
on dividing by t^ — 1 — 2m. 

The existence of this theorem is due to the fact that in z we have the coeffi- 
cients of ^^'' numerically equal and that there are two different sets of equa- 
tions to determine them, since the equation (3e) is not symmetrical with respect 
to f, f~'. We might have obtained another equation, (3e'), in which s, — P 
replace u, P respectively. The sum and difference of (3e), (3e') correspond to 
the rates of change of the areal velocity in the planes of xz, yz. 

Section III. Transformations or the differential equations. 

17. Transformation to avoid the loss of accuracy when r = P {q = 1 or 2). 
The equations (1) are, for terms of order fi, 

Multiply these by s^, u^, 2z^ and sum for all values oi jjl, v except /i = , 
r = , such that 

/A + I' = X. 

The result may be written 

i: \p'(u^s^ + z^z;) - 2Pu^ ■ Ps^ - 2Pz^ ■ Pz^ + 2m{s^Pu^ - u^PsJ) 

in symbolic notation.* 
* Thus the series is 

[("-"»KSl+^("-''»H'^-"»)[Sl+-+^^[31+-L 

where we put u = Uq, s = Sq, z = in the partial derivatives after they have been formed. 
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Next, considering terms of order \, equation (25) may be written 

^^^ =t -n P7] '.+ c- [i>- {I)'/)] : , 

in which the known and unknown terms have been separated. 

Divide the former of these equations by 2 and subtract from the latter so as 
to eliminate the terms in d^J". Also transform the partial derivatives of the 
first line of (8) by means of (6) with /* = . The result may be written 

Dm„ • Z>s, + Z>s„ • Bu^ - {u^D\ + s^I)\) + 2 m (u^Ds^ - s,Z>«.„) 
(9) + £ { - ID\%^.+ ^.^.) + 2(Z>w, • Bs.+Bz^ • Bz,) + m{u^Bs-s^Bu;)} 



n=i 



/-2 



Also, equation (2d) may be written 

+ 2 m (s„ Z>M;, — M, Bs^ + s^ Bu^ - u^ Bs^) 

(10) + glZ)^' (u^s^ + 2^0 - (i)M, • Bs^ + Bz^ ■ Bz;) + 2 m (s^Z)w, - u^Bs;)} 

= i: [(2^ + 1)/Ja + o- iB-^(B'f)y,. 

The terms of principal importance for the determination of X, X' in (9), (10) 
^"■^ _t(X-X')-2(1 +m)(X+X'), 

t(1 + 2m)(X - X') + (t^ + 1 + 4m)(X + X'). 

Now t'=Pj^ when t=P^. Hence, multiplying (9) by (l + 3m), (10) 
by 2 and adding, we obtain an equation for the determination of t (X — X') 
which, when divided by t , i. e., subjected to the operation Z)"' , will give X — X' . 
The equation which results from these processes is, after one or two transforma- 
tions, 

2t(w^s„ + 5;,w„) - (u^Bs^ -f s^Bu^){l - 3m) + 4m(M^s„ - s^u^) 

-6{m/T){u,{B's, + Bs, -mBs^) + s,{B'u,-Bu^ + mBu^)} 

(11) + 3 g 1 1 (1 - m) T (u^ s, + z^ z,) + (m/r) { 2 {Bu^ • Bs, + Bz^ ■ Bz^) 

+ {\-m){s^Bu^-u^Bs,)]^ 

= ^ t [(. + i)/Ja + ^^t f^) P'/]: - '-^--"^ [^'/];- 



174 E. W. BEOWN : ON THE [April 

For inequalities of the first two orders in which t = _Pj (there are none for 
T = I*^) , this equation serves to prevent any loss of accuracy. For in the left 
hand member, t only appears as a divisor in the combination m/r. The first 
term of the right hand member has the combination mYr, the second term is 
absent, the third term has the combination m^r^. 

When T = Pji til® only terms in the left hand member which give rise to a 
loss of accuracy are those multiplied by m/r. For the first of these it should 
be noted that D^s^^ + Z^Sp, D^u^ — Du^ are divisible by m^. 

18. Transformation to ovoid the loss of accuracy when =h t = c + P^ and 
the characteristic contains k^. 
We have the equations* 

(i>^ + 2mZ>).,= [^/]^, (Z>^-2mZ)).,= [|]^, 
(i>^ + 2mZ)).,= [|(]^, iD^-2..D)s^.= Y£\, 

in which b is a characteristic of order X — 2 and not containing k. The right 
hand members of these equations may be expanded in powers of 

U — U^, S — Sj , z, 

the first pair does not contain z , the second pair contains z only in the second 
power, and the third pair only in the first power. The right hand members 
may be therefore written 

"*><' [ dsdu J„ + ^"^ L ^s' Jo ' L Ssdz' J„ ' 

"' L ^«' Jo "' L Ssdu J, + 2 «>= L dudz' J„ ' 

1, m 1. \m +uu r^-^j +^zs \^iA +A'A 



where the functions [J^, []^i,do not contain u^, s^, a^^, and are at least of the sec- 
ond degree with reference to u — u^, s — s^,z or contain the factor m^ . If, then, 

* In these equations it is more convenient to specify the characteristic than its order ; both 
methods are used in this article. 
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we multiply the six equations by s^s, u^,, — s^, —u^, — s^^, a^ respectively, all 
the terms but those arising from []^, [jj^^ will vanish and we obtain 

where Q^ , Qj *^® power series in m — w„ , s — s^, z containing only even powers 
of a and not containing u^, s^, z^^, the first terms being of degrees 3, 2 respec- 
tively. 

For the terms of characteristic bk^ which have the small divisor ± t — c = ^ 
we have to consider the expression (5)'. The portion of this containing m^,, 

+ 2 (., D\, + z,, I)\ + m^ ■ Dz^, ) . 
Adding twice the left-hand member of (12) to this, it becomes 

-2)K.?-)-i)(sj)-Z)Kr->)-i>(vO 

(13) +2D (s^2 Z)m^ + M^2 Ds^ — s^ Du^i — u^ Ds^, ) 

+ 2mZ>(s,,M, - «,,sj + ^z^B'z^^ + ^^\J^\, - 2[^3 + »»' §Jkn>- 

Since the coefficients of f *'^ in this expression are equal, it is sufficient to con- 
sider the coefficients of f '^ . 

First, suppose that t — 2g does not contain g . We have 

«k^r' = (^)o?^^ + (^'')o?-'^ + (i'io), 

where (^^)„, (k' \ are known to m^" (Art. 22), and suppose that 
wJ-i = /3?--F/3'r--F(l,10), 

in which o- -|- 2g = t. Substituting in (13), we obtain for the coefficient of 

/3(^'')„^, 

- 2crg -H 2(<7 + 2g)(cr - 2g -f 2m -f 2) = 2<7(t - g) + (2, 11) = {q, 11), 

since (Art. 20) g = 1 -f m -|- (2, 11), t — g= t — c — (g — c) = (y, 11), 
when T - c = P^ . The coefficient of {¥ X /S' f'' is not small. But {k\ =(2,11). 
For, if we put t = 2g = 2 (1 -|- m) -f- (2, 11), \ = (F)„, \' = {k'^X in (3a) this 
expression becomes 

_(4-4m)-(P)„ + (F^)„P,. 

All the other terms in (2a) either contain the factor m^ explicitly or are formed 
from products (1, 10) (2, 10) and are, therefore, at least equal to (2, 11). 

Trans. Am. Math. Soc. 19 
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Also (k''')^= (0, 10) and, therefore, {¥ )qP^ = (2, 11); as there is no small 
divisor we obtain (P)„ = (2, 11). 
Again, we have 

^,T/-1 = A;„(?« -?-«) + (1,10), 

where ^^ = (0 , 10) (Art. 31) , and suppose that 



Then the term involving yJo^^'^ in (13) has the coefficient t— g=(g',ll) and 
all others terms arising from z^^ are equal to (1, 10) (1, 10) = (2, 11) at the 
worst. 

Collecting these results we see that all the terms of (13) except those in ^s 
are either multiplied by m^, or are obtained to m" or consist of the expressions 
/8(2, 11), /S'(2, 11) — the result which this investigation was undertaken to 
demonstrate. 

Second, suppose that t does not contain 2g . Then 

t,^,r^ = (M')„r + (i,io), 

where (M')„ = (0, 10). "We use the untransformed expression (5'). The first 
two terms then contain the factor (2 , 10) at least and it is only necessary to ex- 
amine (Z>a)' - i)' (2') . 

Section IV. The degrees of accuracy of the coefficients. 

19. The intermediate orbit is the variational curve. It will be assumed that 
the coefficients of the periodic terms in this orbit are obtained to a given degree 
of accuracy (Art. 4) equal to 10, and we shall examine the degrees of accuracy 
to which the terms of the first, second, • • • orders may be found. I use a special 
number instead of a general one ; nothing is lost in generality by this, for if we 
require the degrees of accuracy for any other starting number p , all that is nec- 
essary is to add ^ — 10 to the degrees of accuracy throughout. The special 
number enables me to set down the results more briefly and clearly, and, more- 
over, by Art. 5, they are ready for immediate practical application. 

The argument refers throughout to the power of f in wf~', s 1/ — 1, unless 
otherwise stated. 

It will frequently occur that coefficients may be found more accurately and 
may begin with a higher power of m than the statement shows. In all cases, 
the statement refers to the lowest degree of accuracy and the lowest power of 
m which it is necessary to consider; that is, the expression {p, q) may include 
the cases (p + p', q -\- q) where p' = , §'' = , except in the cases of di- 
visors when p' is always zero and 5' -f §'' is either 11 or 00 . See Art. 7 (v). 

The results are all collected in the table of Art. 31 ; the discussion will be 
more easily followed by referring to this table rather than to the text. 



20. Here u ■= w^ and 
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Terms of Zero Order. 

«o^"" = Zi«i?'' (i = 0, ±1,±2, •••). 

As a^ is the unit of length, we shall put a^, = 1 . Then 

a„=(0,cx,), a. = (2i,10)=(2,10) (i + 0). 

The principal parts of c , g can be found to m" (Art. 4) , so that 
e = l + m + (2,ll), g = l+m+[(2,ll). 

Also the terms containing m" , m' , m^ never all disappear at the same time from 
i + ■i'c + i"g + I'va. , i. e., from any divisor, to the degree of approximation con- 
sidered here (Art. 3). 

Terms of the First Order. 

21. There are terms for which t=P,t— c=P, where q =1, but none 
for which q =2. For such terms the expressions (4) , (5)' respectively reduce to 

— Z>Mj • Ds^ — 2>Sj • 2>w„ + J (mj Bs^ — Sj Bu^ + u^ Ds^ — s„ Du^ ) , 

-D{u,l;-')-D{s,K)-D{u^^-^)-D{s^i;). 

Since 2>(W(|?~'), Z>(S(,f) are equal to (2, 10), these may be written 

(4)i lD{u,s,-s,u^) + {2,l(i)Du, + {2,l(i)Ds, (r=--p,), 

(5X (2, 10)i?(«,?->) + (2, 10)i?(s,?) (r-o=A). 

For the characteristic e' , the coefficients of f ^^^ have the divisor m and there- 
fore all coefficients can be obtained to m'. But (4)j, on division by m, gives, 
for terms other that the principal ones, (1 , 9) (2 , 10) h- m = (2 , 10) , and all 
the other coefficients can be obtained to m^". Noting (3a) , (2cZ) we find that 

7?„,< = (i, 10); ^„ + 7?;, j/.,^: = (2,io) (i+i). 

For the characteristic a, the coefficients of f ^'^ have the divisor 

l-c=(l,ll). 

Treating these coefficients in a similar manner with (5)^, (3a), (2c?) we obtain 

ai>^=(l,10); 3aj^-f a_.^, a. = (2, 10) (i+±i). 

For the characteristic e , there is more difficulty in avoiding loss of accuracy. 
We put 
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and use equations (2a), (2y). One of the pair e^, e^ is arbitrary; put* 

«o - < = 1 • 

The equation will then give e_j, ej = (1, 9), e_j + 3eJ =(2, 10), and all the 
other coefSficients with a degree of accuracy 10 . 

But adding equations (25), (2cZ), we find, for the terms giving u^, 

'o 

from the values already found. Since k/V^ = 1 + 2m + (2 , 10) , we have for 
the coefficient of ^°~^ in this equation, 

((c - 2)^ _ 1 - 2m + (2 , 10)} {e_,. + < + (2 , 10)} = (2 , 11) , 
or, since e_i + e; = (1, 9), (c - 2)^ - 1 - 2m = (1, 11), 

(l,ll)(e_, + e;) + (l,9)(2,10) = (2,ll). 

Hence e._j + ej = (2 , 11) -h (1 , 11) = (1 , 10) . This, with e_j + 3e; = (2 , 10) , 
gives the following results : 

e„,< = (0,ll); e_j,e; = (l,10); 

3e„ + e;, e_j + Sej, e! = (2, 10) (i + o, -1). 

The coefficients of z^ are given in the table. 

Terms of the Second Order. 

22. For the terms for which t= P^ we can use equation (11) which deter- 
mines all such coefficients without loss of accuracy, since the second term of the 
right-hand member does not exist. 

There are no terms for which t = P^. 

The terms for which t — G. = P^orP^ can only occur with the characteristics 
ee', e'a. 

The terms with characteristic ee' are formed of the products of terms with 
characteristics 0, ee' and e, e'; they all contain m as a factor and the small di- 
visors are equal to P^ = (1 , 11). Now the characteristic e' only occurs in (5') 
in the forms D(u^,^~^), D{s^,^), and these are equal to m(l, 10) = (2, 11) 
for the coefficients of f *"". When the multiplication is made with the terms of 
characteristic e, we obtain for the coefficients of ^*°*"»' 



*T. M. M., pt. I, p. 71. 
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(2,ll)(0,10) + (2,10)(l,10) = (2,n); 

this, divided by f^ , gives (1 , 10) . The terms formed by combining the terms of 
characteristics , ee' are at the worst equal to (2 , 10) (1 , 10) = (3 , 11), which, 
divided by (1, 11), gives (2, 10). Hence all are obtained accurately to m*". 

The terms with characteristic e'a are treated in the same way. The terms with 
arguments ± (1 — m) which give a divisor Pj contain the factor m, those with 
arguments ± (1 + m) which give a divisor P^ do not contain this factor. For 
those formed of the products from e' and a we obtain (2 , 10) (1 , 10) = (3 , 11) . 
Considering the combinations , e'a, we find that for arguments ± (1 — m) ,. 
(5)' is equal to (0, 9) (2, 10), which, divided by Pj, gives (1, 9), and that for 
arguments d= (1 + m) , (5)' is equal to (1 , 9) (2 , 10) , which, divided by P^, gives 
(1 , 9) . The remaining coefficients which have no small divisors can be found 
to m^°. Also 

3 (»?«)>, + ('?'«)-K. = (2 , 10) , (va^y, + 3 (v'a),^ = (2 , 10) . 

It is necessary for the third order terms to obtain {v^)i, (»/' )_i to m". That 
this is possible is seen by substituting t = ± 2(1 + m) successively in (Ba); 
the terms of principal importance becomes — 4(1 + m)X', — 4(1 + m)X. All 
the other terms in (2a) are of the form (2, 10) (1, 10) = (3, 11) or are fac- 
tored by m or m^ and there is no small divisor. 

All terms in z are obtained to m'" . 

Thus all inequalities of the second order are obtained to m'" except those 
with characteristic e'a and arguments =b 1 db m ; these are obtained to m'. 

Terms of the Third Order. 

23. It will be shown that all but two coefficients can be found to m'. 

Consider the equations (3 a) , (3/) for coefficients which have no small divisors. 
We shall show presently that in the worst case of small divisors the coefficients 
can be obtained to m* . In the equations for coefficients with no small divisors, 
the coefficients with small divisors are multiplied by (2 , 10) at least. Hence 
no loss of accuracy arises from this source. The only terms of lower orders 
which are not obtained to m'" are the coefficients of f *i=^™ in the characteristic 
e'a. As third order terms are formed in these equations from products of first 
and second order terms, these products in this particular case will be 

(p , 10) (0 , 9) = (1 , 10) unless j5 = , 

and p = only for the coefficients of f '*°. Hence the only coefficients which 
are not obtained to m^" are those of characteristic ee'a and arguments 
± c ± 1 ± m ; of these, the arguments ± (c + 1 d= m) do not give rise to small 
divisors. 
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24. For the case t = P^, the equation (11) is sufficient to avoid all loss of 
accuracy, except for the second term in the right-hand member. But this, being 
of the third degree with respect to powers and products of u — u^, s — s^, z 
must be formed of first order terms only. These are all equal to (1 , 10) at 
least except for the characteristics e*, ek^, which do not give rise to small 
divisors t= P^. On division by t = P^ the results are in any case true to m' . 

The only case in which t = P^ is that of the coefficients of f =^(•'-1-™) in the 
characteristic ee'a. The only term in which we appear to lose more than one 
degree of accuracy is the second term in the right-hand member of (11). But 
this term is formed from the characteristics e, e', a, that is, it is of the order 

(0,10)(1,10)(1,10)-(2,11) = (1,9), 

at the worst. The first term involving m/r in the right-hand member is, at the 
worst, equal to 

m{(0,9)(2,10) + m(0,9)}-(2,ll) = (0,9). 

The last term of the right-hand member involves a loss of one power of m but 
it does not contain terms of characteristic e'a and therefore can be formed to 
m' ; it contains a term involving m~' . This term is cancelled by a similar term 
in the second term of the left-hand member factored by m/r, but as this term 
contains terms from the characteristic e'a it can only be found to m*.* 

25. When t = c + P , the only cases in which we might not be able to obtain 
the terms correctly to m' occur in the characteristics e' a, ek'^, ee' . 

For the only case in which the terms of the first and second orders are not 
known to m'" occur in the characteristic e'a; but the expression (5)' shows that 
these exceptional terms occur only multiplied by terms of characteristic e', that 
is, they are of the form 

(2,10)(0,9)-(1,11) = (1,9), 

since Z)(m^,^-'), -Z>(s^,t) are equal to (2, 10) and g- = 1 for these terms. The 
products of terms of characteristics e'a with those of characteristics e', a do not 
give arguments t ^ c -J- P . 

Also the only terms which contain a divisor P^ occur in the characteristics 
ek^, ee'^, and the terms of lower orders from which these are formed are known 
to m^" . Hence the terms for which q=2 can in any case be obtained to m' , as 
the expression (5)' shows; suppose these (the worst case) are equal to (0, 8)- 
Then for the terms in which 5' = 1 , the only case in which we could not obtain 



* The terms involving m-^ could doubtless be transformed so as to cancel one another before 
development. But I have not thought it worth while to attempt this since this characteristic is 
very small compared to most of those of the third order. It is to be noted that D'f contains the 
factor m^, since t occurs only in /,, in the form f '™ where i* = 2. 
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the results correctly to m' is that which arises in (5') from the combination of 
the terms in ek^, ee' which are equal to (0, 8) with those of order zero. But 
the latter are equal to (2, 10). Hence the worst case in (5)' is (0, 8) (2, 10). 
As this, when divided by (1, 11), gives (1, 9), we are able to obtain the terms 
in ek^, ee'^ for which t = c + Pj to m^. 

The arguments of the terms of characteristic ek^ which are equal to ± c + P^ 
are ± (c — 2g). The combination of characteristics in (5)' are with ek^, e 
with k^, k with ke. The worst case of with ek^ is, since the terms of argu- 
ments ± (c — 2g) do not enter, 

(1,9)(2, 10)^(2, 11) = (1,9); 

they are therefore correct to m' . We use Art. 18 for the other combinations, 
putting b = e . Then 

/3,^' = (0,10), [$3]k = 0. 

since Q^ does not contain terms of characteristics e, ek^, ek, and the terms with 
characteristics k^, k can only arise in combination with e, ke respectively. 
AU the other conditions of Art. 18 are fulfilled and therefore the results can 
be obtained to m*. 

Finally the arguments of the terms in e e' which are equal to ± c + P^ are 
± (2 — c + 2m). The combinations of characteristics in (5') are e', ee' and 
e, e' . The former gives (1, 10) (1, 10) = (2, 11). The only cases in the 
latter which are not also equal to the same, are the combinations of e^, e^ with 
{rf)^, (rj' )_j. But we have obtained (Art. 22) the latter pair to m" and there- 
fore the combinations are (0, 10) (2, 11)= (2, 11). Hence, on division by 
(2, 11) all terms can be obtained to m'. 

26. As there is no loss of accuracy in the equation for z (Art. 15) the only 
terms in z which cannot be obtained to m'" are those arising from e'a and there- 
fore contained in the characteristic ke'a , with arguments ± g ± (1 ± m) ; these 
are equal to (0 , 9) . 

Let* /i (fK+l+°' _ ^-8-l-») + yri' (^K-l-m _ ^-g+l+m-j 

be the terms with argument d= g ± (1 + m) in si/ — 1 , and in equation (2e) 
separate out the terms which are not true to m^". We find 

D {%Dz^,,^ — z^,,^Du^ + u,,^Dz^ - z^Du,,^) = (1,10). 

Equate to zero the coefficients of ^8+i+'"+i , ^g-i-m+i ^^^j isolate the terms which 
are not equal to (1 , 10) . We find 

* So that, according to the general scheme of notation, 
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S{,.-{r,a).^1c,} = {l, 10), {-,.'+ {v'aUk^} = (1, 10). 
Whence 

3;a + /.'= k^{S {rja)y, + (^'a)_;,} + (1 , 10) 
(«) 

= (0,10)(2, 10) + (1,10) = (1,10), 

by the result at the end of Art. 22. 

Let /Xj, fju[ be the coefficients for the arguments ± g ± (1 — m) we then find 
in an exactly similar manner, since 3 (j)'a)y^ + (''?^)->^ = (2 , 10) (Art. 22), 

(a') 3^, + /.; = (1,10). 

The results (a) , (a' ) will be needed later. 

Terms of the Fourth Order. 

27. It will be sufficient, in order to prove the theorem of Art. 4, to show that 
these terms can be obtained accurately to m*. As all coefficients of previous 
orders are known to m' except those of arguments ± (c — 1 — m) in the terms 
of characteristic ee'a, it is only necessary to consider terms which depend on 
these two coefficients and terms which have small divisors equal to P^. 

The characteristic ee'a only occurs in combination with one of the character- 
istics e, e' , a. The small divisors in the combinations ee'a, e and ee'a, a are 
equal to P^ and the terms with characteristics e' , a are equal to (1 , 10) . We 
thus have, for these two combinations, 

(0, 8)(1, 10) ^(1, 11) = (0, 8), m==(0, 8) -(1, 11)=== (0, 8). 

The small divisors occurring in the combination ee'a, e arise from t = c + P . 
Referring to (5)' , we see that such terms are 

which give at the worst 

(c_l-m)(0,8)(0,10) = (2,10), (1,9 (1, 9) = (2, 10), 

and these, divided by (2 , 11) give (0 , 8) . 

It is therefore only a question of examining terms for which t = P^ or 
r = c + P^, leaving aside the parts containing the characteristic ee'a. In all 
cases these terms can be obtained to m'^ . Hence the contributions of terms of 
orders , 4 will be at the worst 

(2, 10)(0, 7) - (1, 11) = (1, 8), (2, 10)(0, 8) - (2, 11) = (0,8), 

so that these combinations need not be further considered. 
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28. For the terms in which t= I*^, the equation (11) is used. We have just 
seen that the terms with characteristics ee'a will not occur, and therefore the 
only portion which can produce results not correct to m' will be the second 
term of the right-hand member. But this being of the third degree in u — u^, 
s — Sfl, 2 can only contain terms of the first and second orders ; with the excep- 
tion of terms of characteristic e'a, these are all known to m'°. But these ex- 
cepted terms can only occur with terms from two of the three characteristics 
e,, e', a, so that at the worst they occur in the combination 

(0,9)(0,10)(1,10)^(2,11) = (-1,8). 
The terms in m~' disappear and all the results are known to m". 

29. The only terms in which r = c + P^ occur in the characteristics 

e'' a, e'a'\ e'k^a, e^e'a. 

The combinations of orders 1,3, excluding e, ee'a, are all at least equal to 
(1 , 10) (0 , 9) = (1 , 10) so that on division by (2 , 11) they give the results cor- 
rect to m". The only combinations of orders 2 , 2 which do not give results 
correct to m' are e'a, k" and e'a, e^, for all the others are at least equal to 

(0, 10)(0,10) = (0,10). 

For the combinations e'a, k- and k, ke'a, with arguments ± (2g — 1 — m), 
put b = e'a in Art. 18. Then 

/3, /3' = (0,9), Q^ = (0,10), 

for Q^ can only contain terms of the first two orders from amongst which those 
of characteristic e'a are excluded. All the other conditions are fulfilled and 
therefore the results are correct to m^. 

For the same combination with arguments 1 + m , we see from the result at 
the end of Art. 18 that it is only necessary to consider 

for the combinations of arguments — g, g + 1 -f m and g, — g + 1 -f m since 
these are the only ones which will not be equal to (1 , 10) . Using for the 
moment, the notations of Art. 26, these terms give 

;g (g + 1 + m) -f (1 -f m)-} fik„ + {si- g + 1 -f m) + (1 -f m)-} fi'k^ 

which is equal to 

(3m -f m')/.\, + (1, 11)(0, 9)(0, 10) = (1, 10)(0, 10) -f (1, 10) = (1, 10), 
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since (Art. 26) Sfi + fi' = {1, 10). On division by (2, 11) the result will be 
correct to m^ 

For the combination e^, e'a and arguments ± (2c — 1 — m), the only terms 
in (6)' which will not be true to m*" will be formed from 

Substituting in (5)' and neglecting a factor (0, 11), we find that the portion of 
the coefficient of f ^o-i-m ^i^Jcji arises from these terras is 

{^\My. + (^")o(V«)-^ = mo - 3 (e")J (^a)^ + (6'^), {3 (r;a)^ + (Va)_^}. 

It is easy to show* that {e^\— 3(e' )^ = (1, 10), and therefore using the re- 
sult at the end of Art. 22 , the expression is equal to 

(1,10)(0,9) + (0,10)(2,10) = (1,10). 

On division by (2, 11) the result is true to m'. 



30. For terms of orders higher than the fourth a similar treatment can be 
made. If we had obtained all the second order terms correctly to m^" there 
would have been no difficulty in obtaining all the third order terms in u to m' 
and those in z to m^" and then no special investigations would have been re- 
quired to obtain the fourth order terms in u to m' and in « to m' . In fact, the 
transformations of Section III enable us, in the worst cases of small divisors 
(t = J-'j , ± T = c + P^ ) , to throw the loss of two powers of m on the order 
next but one below that under consideration. Suppose that we know the terms 
of order X (X > 3) in u correctly to m''"'^ and it be required to show that the 
terms of order X -f- 1 may be obtained correctly to m'~'^~^ . In the case of these 
divisors the loss of two powers of m is thrown on to the terms of order X — 1 
and lower; these being known to m*~'^+^, the results will be true to m'"'^"*. 
As the third and fourth order terms follow this rule, terms of all the following 
orders will do so. 

The terms of order X + 1 in » can be obtained to the same degree of accuracy 
as the terms of order X iu m , that is to the order m*""*^ . 

The terms of the fourth order being all known to m^ and as we require those 
of the fifth order only to m'' and not to m*, the great majority of the difficul- 
ties arising from the terms of the second order in e'a will not occur after the 
fourth order. 

No general transformation has been given for terms in which ± t = c -f P^ • 
These terms must contain some power of e, e', a in their characteristics; the 



* The expression is zero in the purely elliptic formulie (m = 0) and each coefficient is known 
to ni'". 
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only cases where difficulties will occur will be when they contain e. Such dif- 
ficulties can, however, be solved by the method used in Art. 25 for the terms of 
arguments d= (2 — c + 2m). 

31. The following table shows the results of this section in a form which is 



convenient for reference, 
readg + P . 



In the case of the coefficients in zv^ — 1, for c + P 



Ord. 


r^P„ ±r + c+P, 


r = P, 


±r = c+A 





oo = !.«*= (2, 10) 


none 

7„,'?;=(i,io) 


none 


1 


'oi ^0' *^o=(0, 10) 
otherooeff. = (2, 10) 


e-i,e[, k-i=(\, 10) 
e_i + 4 = {2, 10) 


2 


(^Mo. (^'\, («')o = (0, 10) 

(fc'')„, (tt')o=(0, 10) 

(£2)„-3(£'')o=(0, 10) 

{A;%(r,»)-i, ('/')-!= (2, 11) 

other coeff. =(1, 10) 


all coeff. =(1, 10) 


(?/a)i^, (v'a)-y, = (0,9) 
{va)-y,, {v'a)y,= (l,9) 
3{va)^^+(v'a)-y,= {2,10) 
{1a)-}i + {3v'a)^y,= {2,10) 
other coeff. = ( 1 , 10) 


3 


all coeft. = (0, 9) 

3{k7ia)y,+ {kr/a)^y^ = {l, 10) 
Hkv'a)x+ikva)-y, = {I, 10) 


Other coeff. = (0, 9) 


all coeff. = (0, 9) 


4 


all coeff. =(0, 8) 


all coeff. = (0, 8) 


all coeff. =(0, 8) 



c— 1— ra=(2, 11), 



•1— m = 



(2,11). 



32. The results of this section once proved, the practical difficulty of carrying 
them out is not great. Suppose that it is a question of inequalities such that 
± T = c + P. Then 3X + X' can be and would be determined without loss of 
accuracy and without any increase in the number of calculations, although there 
may be a loss in X, X' . In the after work we simply treat the result for X — X' 
as though it were accurate (if a numerical development is being made, an arbi- 
trary figure is added to the last place of decimals found) ; it may be the func- 
tion 3X -f X' which afterwards really arises. If it be so, the future results 
will be correct to the orders stated above. 

When this method cannot be followed, it is usually best to find the coefficients 
as accurately as possible by the ordinary methods and then treat the special 
coefficients in which there has been a loss of accuracy by the special methods. 
The case of €_, , e^ (Art. 21) is an example of this. 

Haveeford College, Pa. 



